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Abstract. We describe random loop models and their relations to a family of quantum 
spin systems on finite graphs. Spin correlations are given by loop correlations. Decay of 
correlations is proved in 2D-likc graphs, and occurrence of macroscopic loops is proved in 
the cubic lattice in dimensions 3 and higher. A phase transition with magnetic long-range 
order is proved for SU(2)-invariant models of spin-1. 
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We study models of random loops that are closely related to quantum lattice systems. 
There are many reasons to consider them. First, the probabilistic setting is mathematically 
elegant, it is fairly simple to depict, and the presence of spatial correlations makes it very 
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interesting. Second, quantum spin systems play a major role in our understanding of the 
electronic properties of condensed matter. And third, the relations between these models 
are fascinating. They hold in spite of their mathematical differences, and several properties 
happen to be visible in one representation but not in the other. 

The models of random loops described in this article find their origin in the work of Toth 
on the quantum Heisenberg ferromagnet |54j , and in the work of Aizenman and Nachtergaele 
on the Heisenberg antiferromagnet [3]. The present extension allows to describe the spin-| 
quantum XY model and certain SU(2)-invariant spin-1 spin models. A major advantage 
of these representations is that correlations between quantum spins are given in terms of 
properties of the loops. 

Quantum Hamiltonians are also related to the generators of evolution of classical inter- 
acting particles. This is an active research area and some methods have been borrowed 
from quantum systems, such as the Bethe ansatz. We refer to [531 [23 EH] for recent 
discussions. Further connections should be possible and the probabilistic representations of 
quantum lattice systems at equilibrium could help to bridge the gap. 

Many important results have been obtained for quantum spin systems that are relevant 
for models of random loops. One such result is the Mermin- Wagner theorem about the 
absence of continuous symmetry breaking in dimensions 1 and 2 [ISJ US EH \EH HH 123 031 
[351 [23 [HUES]- We provide a simple proof; it uses operator theory and it is restricted to those 
parameters for which there exists a correspondence with quantum models. An interesting 
challenge is to find a probabilistic proof that holds more generally. 

Another important result in statistical mechanics is the proof of existence of a phase 
transition with continuous symmetry breaking and long-range order in dimensions 3 and 
higher. The first positive result is due to Frohlich, Simon, and Spencer, for the classical 
Heisenberg model |26j . The extension to quantum models was achieved by Dyson, Lieb, 
and Simon |19j , and it is of direct relevance here as it implies the occurrence of macroscopic 
loops. They used the method of reflection positivity and infrared bounds, that has been 
subsequently discussed and extended in [50l E21 [231 EH [42j El 0] . 

The main goal of this article is to show that this method can also be applied to the models 
of random loops. As a result, we prove the occurrence of macroscopic loops when the spin 
parameter S is small and the dimension of the lattice is large enough. In the case S = we 
recover the results of [HI [41] . The present approach becomes useful, from the perspective 
of quantum spin systems, in the case S = 1 since it allows to uncover the existence of long- 
range order that is neither ferromagnetic nor antiferromagnetic. Besides, the representation 
of random loops suggests that correlations of the form (S^Sy) decay exponentially fast with 
respect to \\x — y\\. 

The model of random loops is introduced in Section [2j We discuss the family of quantum 
systems in Section[3]where we state and prove the connections to random loops. Correlations 



for SU(2)-invariant systems are particularly interesting, see Theorem 3.5 The Mermin 



Wagner theorem is discussed in Section [4j We prove it using a simple inequality that, in 
essence, estimates the difference of expectations with respect to two Gibbs states in term of 
their relative entropy. The occurrence of macroscopic loops is describe d in Sec tion s [5] and 



tWe discuss explicit models with spin | and 1 in Section u\ 
e new results for quantum systems obtained in this article. 



7.2 



and 



7.3 



contain 



Theorems 

the new results for quantum systems obtained in this article. We conclude with a heuristic 
discussion about the joint distribution of the lengths of macroscopic loops, and of the nature 
of pure Gibbs states of the quantum models. 
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2. Probabilistic models 



We start by introducing the model of random loops in Section 2.1 then the model of 



space-time spin configurations in Section |2.2| The former is more elegant and it contains 
all the relevant events for the description of quantum correlations. The latter contains more 
information and it will be used in Section [5] that discusses the occurrence of macroscopic 
loops in rectangular boxes. The model of random loops is a marginal of the second model. 



2.1. Model of random loops. Let (A, £) denote a finite graph, with A the set of vertices, 
and £ the set of edges. An edge is always between distinct vertices, and two vertices are 
connected by at most one edge. The most relevant example is the box A = {1, . . . , L} d in 
Z d , with £ the set of nearest-neighbors, £ = {{x, y} C A: \\x — y\\i — 1}, but we consider 
more general graphs. Let j3 > be a parameter. We consider a Poisson point process on 
£ x [0, j3] with two distinct events, the "crosses" and the "double bars" . Let Wk denote the 
collection of subsets of £ x [0, 0\ with cardinality k. The set ft of the realizations of our 
Poisson point process is then 

n = \J w kl x \J w k2 . (2.D 

fci>0 fe 2 >0 

Let u G [0, 1] be another parameter. We let p u denote the probability measure on 17 that 
corresponds to Poisson point processes of intensity u for the crosses and 1 — u for the double 
bars. 

To a given realization ui G of the Poisson point process corresponds a set of loops, 
denoted C(oS). The notion of loops is best understood by looking at pictures, see Fig. [T] It 
is a necessary pain to go through mathematically precise definitions, though. So a loop of 
length £ is a closed trajectory 7 : [0, f3£] pcl - —> A x [0, /3] por such that 

• 7( s ) 7^ 7( s ') if s i s ' are distinct points of continuity of 7. 

• 7(5) is piecewise differentiable, and 7'(s) = ±1 at points of differentiability. 

• If s is a point of non-differentiability, then {7(5— ), 7(s+)} G £ x [0,/3]. 

We identify loops that have the same support but different parametrizations, i.e., the loops 
described by 7(s + t) and j(—s) are considered to be identical to 7(5). 








Figure 1. Graphs and realizations of Poisson point processes, and their 
loops. In both cases, the number of loops is = 2. 
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Let uj € 0, and (x, t) 6 A x [0, /3], but not in the support of oj. The loop 7 that contains 
(x,t) can be defined by starting with "f(s) = (x,t + s) in a neighborhood of s = 0, and by 
moving "up" until a cross or a double bar is met. If it is a cross, go across and continue 
in the same vertical direction. If it is a bar, go across and continue in the opposite vertical 
direction. The trajectory will eventually return to (x,t) and the loop is complete. Let C(uj) 
be the set of all loops, starting from any point of A x [0,/3] por (not in the support of oj), 
modulo reparamctrization. The number of loops is necessarily finite (it is always less than 
|A| + M). 

Let 9 > be a parameter. We define the partition function as 

y e (M) (/?,A)= / e^d Pu (uj). (2. 2 ) 



12 



The relevant probability measure for the model of random loops is then given by 



e\ C ^d Pu (cj). (2.3) 



Correlations will be given by three events, i.e., subsets of f2: 

• E xy t is the set of all w € il such that (x, 0) and (y, t) belong to the same loop, and 
with identical vertical direction at these points: 



d 
ds 



7(s) = T-70) 



(2.4) 



7(s)=(x,0) ds 

(This definition does not depend on the actual parametrization of 7.) 

E~ yt is the set of all u G il such that (x, 0) and (y, t) belong to the same loop, and 

with opposite vertical directions at these points: 

-lis 



(2.5) 

i(s)=(y,t) 



f(s)=(x,o) ds 

• Ex,y,t — E xyt U E~ y t is the set of all realizations uj € f2 such that (x, 0) and (y, t) 
belong to the same loop. 
We also define the length L x j of the loop that contains (x, t) € A x [0, j3] as the sum of all 
its vertical lines. 

The case u — 1 and 9 = 1 is the "random interchange model" , or "random stirring" (35] , 
and loops are closely related to permutation cycles. The joint distribution of cycle lengths 
was obtained by Schramm for the complete graph [55 following a conjecture of Aldous; see 
also [7] . Infinite loops have been proved to occur on trees [331 IS] ■ Exact formulae for the 
probability of cyclic permutations have been obtained in There is no correspondence 

between random loops with = 1 and quantum systems, and, rather curiously, none of 
the results obtained in the present article (decay of correlations in d = 2 and occurrence of 
macroscopic loops in d > 3) apply in the seemingly simpler situation 9 = 1. 

It is not hard to prove that, when (3 is small, V(E X}Vtt ) decays exponentially fast with 
respect to the distance between x and y. See e.g. Theorem 6.1 of [25] , 

2.2. Space-time spin configurations. Representations with space-time spin configura- 
tions have proved useful in many instances, see e.g. [2H1 HOI 113 HH 1301 |3B] 130] [T7] and 
references therein. Given S € ^N, a space-time spin configuration is a function 

a : A x [0, /3] pcr — > {-S, -S + 1, . . . , S}. (2.6) 

such that cr x t is piecewise constant in t, for any x. See Fig. [2] for an illustration. Let 
E denote the set of such functions with finitely-many discontinuities. We further require 
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that these functions be cadlag so they can be equipped with the Skorokhod topology. Non- 
experts should pay no attention to these technical details, they just help define a space of 
well-behaved functions on S. 



■l +1 
o 



-i 





+1 

-1 



+ 1 
-1 

O 



Figure 2. Illustration for a realization of the process p L , for which there 
exists a space-time spin configuration. 



If u> is a realization of the Poisson process p u described above, we say that a is compatible 
with u if o~ x t is constant on each loop of C(u>). Let denote the set of all such space-time 
spin configurations. Notice that 

|E(w)| = (25 + l) |£(w)l . 
For 8 = 2S + 1 we have the relations 



(2.7) 



r e (u) (/3,A)= I dp u (u) 



1 



and 



lS(S + l)V(E x , v , t ) = 



1 



Y 



(«)/ 



The latter relation follows from the following identity, which holds for all S G |N, 

1 



s 



2S 



1 ^ 

a=-S 



a 2 = ±S(S + l). 



(2.8) 



(2.9) 



(2.10) 



We introduce now a more general setting for measures of space-time spin configurations. 
This will help us to make the connection between quantum spin systems and random loops. 
More importantly, we shall need the flexibility of this setting in Section [5] in order to prove 
the existence of macroscopic loops. 

It is enough for our purpose to consider only discontinuities that involve nearest-neighbor 
vertices. We introduce a Poisson point process on £ x [0, (3] whose objects are "specifications" , 
that specify restrictions on the local configurations. Namely, let {x, y} £ £ and t £ [0, /?]. 
An object is a set A of allowed configurations at x and y, immediately before and after t. 
The configuration at this "transition" is conveniently described as follows 
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For this notation to make sense we need an order on A; in the configuration above, we 
assume that x < y. A set A is a subset of {— S, . . . , S} 4 and it appears with intensity i(A). 
We are now ready for precise definitions. 

Let i : V({-S,...,S} 4 ) — > M.+ denote the corresponding intensities (V denotes the power 
set) and p L the Poisson point process. Given a realization £ of p L , we let £(£) C S denote 
the set of compatible space-time spin configurations, that is, a € £(£) if 

€ A whenever £ contains the set A at (x,y,t). 



• er T * is constant in t otherwise. 



>x,t 

The measure is then given by p L and the counting measure on compatible configurations. It 
is important to remark that different intensities i and if can give the same measure. They 
can be characterized by their "atomic composition" a,, : {—S, . . . , S} 4 — > R+: 



«.(H)= £ **>■ 



(2.11) 



A 3 = 

c d 

Lemma 2.1. We have 



J d Pl (o e F (°) = J d ^) E F ^ 



/or a// -F €E C(E), «/ and on??/ i/ a t = av . 

It is not hard to prove this lemma, e.g., by discretizing the interval [0, 0\ and by comparing 
the sums on both sides. An immediate consequence is that Poisson point processes can be 
composed as follows. 

Lemma 2.2. Let i, t! be intensities V({—S, . . . , S} 4 ) —> K+. Then 

f dp 4 (0 f dpAC) E F l?) = [ d/w(0£*V) 



/or any F G C(S). 



Indeed, this follows from a t + a,/ = Next, we relate the present setting with the 

random loops of the previous subsection. 

Lemma 2.3. Assume that 9 = 2S + 1, and let the intensity i be defined as follows: 



V L ft a J a,6G{-S,...,S}/ Via a J a,b£{-S,...,S}/ 

We set l(A) = otherwise. Then 



l-u. 



and 



|5(5+l)P(^, t )= 1 / dp t (0 E 17 

*6 (P) ^) — 



This follows from Eq. (2.9), once we observe that the sets above describe precisely the 
specifications given by the crosses and the double bars. 
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3. Quantum spin systems 

3.1. Families of quantum spin systems. Let S G |N. The Hilbert space that describes 
the states of the system is the tensor product 

n A = <g)n x , (3.D 

xeA 

where each H x is a copy of C 2S+1 . Let S 1 , S 2 , S 3 denote the usual spin operators on C 2S+1 . 
That is, they are hermitian matrices that satisfy 

[S\S 2 ] =iS 3 , [S 2 ,S 3 ] = iS\ [S 3 ^ 1 ] = iS 2 , (3.2) 

(S 1 ) 2 + (S 2 ) 2 + (S 3 ) 2 = S(S + l)ld. (3.3) 

Recall that each S l has spectrum {— S, -S + 1, . . . , S}. We use the notation S = (S 1 , S 2 , S 3 ) 
and, for neR 3 , 

S s = a-S = a 1 S 1 +a 2 S 2 + a 3 S 3 . (3.4) 
These operators are related to the rotations in R 3 by 

a -iS 3 ab AS 3 cR s b ,„ . 

e o e = o (3.5) 

where Rgb denotes the vector b rotated around a by the angle ||a||. 

Let S % x — S l ® IdA\{ x }. We use Dirac's notation since it it very convenient. In C 25+1 , \a) 
denotes the eigenvector of S 3 with eigenvalue a £ {— S, . . . , S}. In H x (8) H y , \a, b) denotes 
the eigenvector of both S 3 and S 3 with respective eigenvalues a and b. 

We consider the three operators T xy , P xy , Q xy on / H^ x y y (and their extensions on T-La by 
identifying T xy with T xy ® Id A \{x,i/}, etc ---) : 

• T xy is the transposition operator: 

T xy \a, b) = \b,a). (3.6) 

• P xy is the operator 

s 

P X y= (-l) a - b '\a,-a)(b,-b\. (3.7) 

a,b=-S 

Equivalently, the matrix coefficients of P xy are given by 

(a,b\P xy \c,d) = {-l) a ~ c 5 a ^ b 8 c ^ d . (s.s) 

Notice that 2g 1 +1 P xy is the projector onto the spin singlet. 

• Q xy is identical to P xy except for the signs: 

(a, b\Q xy \c, d) = 6 a ,bS c ,d- (3.9) 

The first two operators are invariant under all rotations in M 3 , the last operator is invariant 
under rotations around the second direction of spins, as stated in the following lemma. 

Lemma 3.1. For all oel 3 , we have 
(&) c ' x ' y P xy e ' x ' y T xy . 

(b) e - iS *- is3 yP xy e iS * +is3 y =P xy . 

And for a — ae*2, a <E M; or a = a\e\ + a 3 e 3 such that a\ + a\~ ir 2 , 

(c) e-^-^Qzy =Q xy . 
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Proof. For (a), it is enough to check the following commutation relations: 

[Si + S l y ,T X y] = 0, (3.10) 

i = 1, 2, 3. Commutation with S x +Sy is obvious. It is convenient to introduce 5* = S^iiS 2 . 
It is well-known that 

S ± \a) = y/S{S + l) - a(a±rj\a± 1), (s.n) 
with the understanding that 5' + |5') =0 and 5~| — 5) =0. Then it is not hard to check that 

[S± +S y t ,T xy ]=0. (3.i2) 

It is not easy to check (b) directly. But 2 s+i ^ x v IS tne projector onto the eigenspace of 
(S x + S y ) 2 for the eigenvalue 0. The eigenspace is known to have dimension 1, see e.g. [35], 
and the result follows. Finally, (c) follows from the relation 



and 



e" i7rS » P xy e i7TS y = e^ s y P xy e~ i7lS y , (3 .iaj 



e --(^+S?) e -- s » e i < s 'i +s ^ = e l < . (3 .i4) 



□ 

We consider two distinct families of Hamiltonians, indexed by the parameter u £ [0, 1]: 

(3.15) 



{x,y}££ 

Ha ] = - ( uT *v + u ) p *y) • ( 316 ) 

The first family is convenient for the probabilistic representations, and it contains many 
relevant special cases for S = , : The usual Heisenberg ferromagnet and antiferromagnet 
models, and the XY model. This is explained in Section [7Tlj The second family is physically 
more relevant and the S = 1 will be treated in details in Section 17.21 

Let Z( U \(3,A) and Z^ u \j3, A) denote the corresponding partition functions: 

zW(/3,A)=Tr WA e-^i"\ ( 3 .ir) 
Z»(/3,A) = Tr„ A e-^ U) . C 3.i 8) 

The expectation of the operator A in the Gibbs state with Hamiltonian is 

(A) = . .} , Trie-^" . (3.19) 

We also consider the Schwinger functions; given t £ [0, j3] and two operators A and B on 
Ha, let 

{A;B){t) = z(^^) TlAc ^ t)<} Be " <} • (3 ' 20) 

Notice that (AB) = (A;B)(0). 
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3.2. Random loop representations. The representation of the Gibbs operator e~@ H 
in terms of probabilistic objects, with H a Schrodinger operator, goes back to Feynman's 
approach to the interacting Bose gas. Such representations of lattice systems has allowed 
many authors to prove the occurrence of phase transitions in anisotropic lattice systems 
[28l l40l H5l H8] . Conlon and Solovej used a random walk representation in order to obtain 
estimates on the free energy of the S = \ Heisenberg ferromagnet [16] . Their result was 
improved by Toth using a loop representation [54 , . A similar representation was introduced 
by Aizenman and Nachtergaele for the S = | Heisenberg antifcrromagnet, and more gen- 
erally for interactions of the form P xy [3]. It allows them to relate the quantum spin chain 
(d = 1) to two-dimensional Potts and random cluster models. Nachtergaele has proposed 
extensions for higher spins in [47l [48] . In this section we show that the representations 
of Toth and Aizenman-Nachtergaele can be combined and extended to the families 



and defined in Eqs (|3.15|) and (|3 . 16h . As it turns out, the representation holds for all 



r(«) 



S E in the case of the family but only for S E N in the case of the family 

The first identity between quantum system and loop model concerns the partition func- 
tions. 



Theorem 3.2. For all u E [0, 1], we have 



(25+1) c ^'dp u (uj) = > ~ ' " ' " 1 



ZW(jS,A) for all S e §N, 
Z(")(/3,A) for all SEN. 



Proof. Using Trotter's product formula 

' -in. Tr( H e A 

N 



Tr e"^ " = e^ £ I limTr TT e U^+{^)Q* v -i) 



N 



" " Hm Trf [] [l - § + § (uT xy + (1 - u)Q a 



N—fOD 



{x,y}e£ 



N 



N 



* m J™ E R(^~£ + £(uT xy + (l-u)Q xy )]\a^). 



T (i) cr< N ) i=l 



(3.21) 



The sum is over <jW g {—S, . . . , S} A and we set a( N+1 ' > = The transposition operator 
yields the specification 

f — i 

lb a J a,be{-S,...,S} 

and the operator Q xy yields 

La a J a,be{-S,...,S} 

In the limit N — > oo we obtain the expression of Lemma |2.3| This proves the claim for 

The proof for Z^((3, A) is similar, except for two differences: 

• In order for the product of matrix elements to differ from 0, the spin in the loop 
must change sign when the loop changes its vertical direction (that is, at double 
bars). 

• The matrix elements of double bars are 



(-1) 



<7x j l± 3 <Tx j '*J 
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and they can be factorized with respect to the loops: (— l) (Tx r t j- for the loop coming 
from below, and (— l) <Tx j t i+ for the loop coming from above. If the spin of the loop 
is even, all factors are equal to 1. If the spin is odd, the factors are all equal to — 1 
and their product is 1 because there is an even number of them. 

□ 

The situation with half-integer spins is very different. Some loops receive negative weights, 
such as a loop with two neigboring sites and two transitions, one cross and one double bar 
(see Fig. [3| . The representation therefore involves signed measures (still real) and we lose 
the probabilistic nature. It is not clear whether such representations can be useful. 




1/2 



-1/2 



6 6 



Figure 3. A "bad loop" on two vertices for the representation of the 
family with 
(-!)!-(-!)= -1. 



family with half-integer spin. The factor due to the double bar is 



We turn now to correlation functions. It is remarkable that the spin-spin correlations 
of the quantum models can be expressed in t erms of properties of loops. Recall the events 
E^yj introduced above (with the help of Eq. (2.4|). We state first the results for the family 



H^; the results for H 



rr{u) 



3.5 



i A , unc irouiio im _i_i A are postponed to Theorem 

Theorem 3.3. Consider the Hamiltonian with S G and u G [0, 1]. Correlations 
the spin directions 1 and 3 are given by 

{Sl-Sl)(t) = {Sl-Sl){t) = lS(S + l)V(E x , y , t ). 

Correlations in the spin direction 2 are given by 

{S 2 X - S 2 y )(t) = |S(S + 1) [V(E+ Vtt ) - P(B- Vlt )] ■ 

Proof. Equality of correlations in directions 1 and 3 is clear by symmetry. We have 

1 



(S 3 x ;S 3 y )(t)- z(u)(/3jA) 



£(«)(/?, A) 



/ dp u (uj) E a x,oo- y ,t 



(3.22) 



a— — S 



The result follows from Eq. (2.10). 



For the correlations in the direction 2, we write a similar expansion as in Eq. (3.21) but 



with additional factors (a 



x,0- 



■;0+, 



and {(T y .t-\S x \o-yj+)- These factors force (x, 0) and 
(y, t) to be in the same loop. Now recall that S 2 = ±(S+ - S~) while S 1 = + 5"). If 

ui G E X y t> there is one factor with S + and one factor with S*~, resulting in — i 2 times the 
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same contribution as for S 1 . If u e E~ y t , on the other hand, both factors involve S + or 
both involve S~ , and the contribution is i 2 times that of S 1 . □ 

Macroscopic loops are related to two physical properties of the system, namely sponta- 
neous magnetization and magnetic susceptibility. This is stated in the following theorem. 

Theorem 3.4. 

(a) Macroscopic loops and magnetic susceptibility: 

3 d 2 , „ _„„(«)_ 



xGA 



h=0 



(b) Macroscopic loops and the expectation of the square of the magnetization: There 
exists a constant K (it depends on S but not on (3,u,A,£) such that 

-|^- e (s^)-k^m^)V\mW\<k(J2l ( ^) < ^l^y E <^ 3 >- 

x.yeA xeA ^ ' x.yeA 

The claim (a) states that, for regular domains such as cubes, the infinite volume free 
energy is non analytic when the average length of the loops diverges. In the claim (b), the 
difference between the upper and lower bounds is of order |A| (if the graph has bounded 
degree) and it is much smaller than E(^i( a , )) when loops have infinite average length. 

Proof. Let us introduce the Duhamel two-point function product: 

An) 



1 fP 

[ ' jDuh - ZOO (/?, A) J Q a e 



(3.23) 



Using Duhamel formula, we have 

Tre-^'+^E^A^ _ ^ U )(I3, A) |l + \j3h 2 E (S 3 X , S 3 y ) Duh + O^)]. (3.24, 

x,y£ A 

Notice that the odd powers of h do not contribute because of symmetry (rotation around 
S 2 of angle n) . We have 

E (E i (^° ) ) = E fnE x , y ,t)dt 

x£A x,yeA 

X.J/fcA (3.25) 



x,yeA 

^logTre-^+^E^A^ 



h=0 



f3S{S + 1) dh 2 
This proves (a). 

Let us recall the following inequalities that relate expectations with respect to Gibbs 
states and the Duhamel two-point function: 

±(A,A) Duh <±(A*A + AA*) 

. / (3-26) 

< ^(A,A) Uuh ^([A*,[H^\A}}) + j(A,A) Buh . 
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The first inequality follows from the convexity of the function F(s) = Tr A* e~ sH ^ A e~(@~ s ^ H A 
The second inequality is more involved and is a consequence of the Falk-Bruch inequality, 
that was proposed independently in [5D] and [H] . Notice that the expectation of the double 
commutator is always nonnegative because it is equal to ([H^\A], [H^ \ A]) Du h- Using the 
inequality above with A — we immediately get the upper bound in (b). We need 

to deal with the double commutator for the lower bound. It can be calculated, and we will 
actually need it in Section [5j Here, it is enough to notice that 



+ Sy 



[T X y , S x 



S 3 y}]=0, 



(3.27) 



so that 



{x,y}e£ 



(3.28) 



x,y 



where Q x 3 is defined below in Eq. (5.39). We only need here that Q xy be bounded, which 
is obvious. Inserting in the second inequality in (3.261, we have 



]T (S 3 x S 3 y ) < K E(J2 h x , 0) ) VTl - u)\£\ +E{J2 i(x,o) 



(3.29) 



Using E(%2 X L( x )) < (3\A\, the result follows. 



□ 



Let us turn to the family of rotation invariant Hamiltonians H A . 

Theorem 3.5. For the family Hj^ with SeN and u G [0, 1], we have the following relations 
between spin and loop correlations: 

(a) (su s*)(<) - is(s + 1) [mt.y.t) - nK.y.t)] ■ 

(b) ((Sif; (S;)2)(t) - ((^) 2 )((^) 2 ) = ±S(S + 1)(2S - 1)(2S + 3)F(E x , Vft ). 

The claim (a) suggests that correlation such as (S x S y ) decay fast for u 6 (0, 1). This will 
be a striking feature of the model with spin 5 = 1 which is discussed in Section [772] Notice 
that if u = 1, E~ y t has probability zero. And if u — and the graph is bipartite, it is not 
too hard to check that E xy t , resp. E xy t , has probability zero if x and y belong to different 
sublattices, resp. identical sublattices. 

The constant in the right side of (b) is zero for S — \ but it is positive otherwise. 



Proof. The model is rotation invariant, so it is enough to prove it for i — 3. The correlation 
function of (a) is non-zero only if a loop connects (x,0) and (y,t). Spin values are identical 
if the vertical direction in the loop is identical, they are opposite if the direction in the loop 
is opposite. Hence the difference of probabilities of E^ y t and E~ t . 
For (b), let us first observe that 



<(s;) 2 ) = l((sir + {si? + (S 3 x f) = ±S(S + 1). 



(3.30) 
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Next, 



= ^Ay4,., d -<")< 2S+1 > |£MI (5^r£ 

dM „ )(2S + 1) ^»i^L_ £ 



S N 2 

a 2 

a=-5 
S 



(3.31) 



;f5(5+l);r 



a 4 



^ £ a 4 - (15(5 + 1)) 



a=-S 

2 1 



25 , 

a= — b 



We used the identity ( |2.10[ ). The claim follows from the identity 



1 

— ]T a A -{\S(S + l)) 2 ^^SiS + mS-mS + 'i), 



25 , . 

a— — b 



(3.32) 



which is valid for all integer and hall-integer 5. □ 

4. Decay of correlations in 2D-like graphs 

It is well-known that continuous symmetries cannot be broken in two spatial dimensions. 
This was first proved by Mermin and Wagner for the quantum Heisenberg model [45 . The 
original statement is about the absence ol spontaneous magnetization in the lattice Z 2 . 
Many improvements have been made over the years, in particular the prool by Frohlich and 
Pfister that all KMS states are rotation invariant [231 [25] . Here we are especially interested 
in the decay ol the two-point correlation lunction. The first result for the Heisenberg model 
in Z 2 is due to Fisher and Jasnow [H] . Algebraic decay was proved by McBryan and Spencer 
for the classical Heisenberg model, in a short and lucid article that introduces the idea ol 
"complex rotations" [33] . This was extended to quantum systems in [T3l [35] — the bounds 
are given in the form ol an integral ol Bcssel lunction, and algebraic decay could be extracted 
from it. Further work on the ground state ol systems with less than two dimensions can be 
found in [33] . 

Here we re-derive the result ol |21| in a simpler lashion. The present prool is inspired 



by [52 1231 GUI [3S]- Theorem 4.1 is also more general as we do not assume that the graph 
is regular. A consequence is the absence ol macroscopic loops in 2D-like graphs for 8 — 
2, 3, 4, ... . As the prool relies on the continuous symmetries that are present in the quantum 
setting, it does not seem possible to extend it to other values ol 9. 

Let d{x,y) denote the graph distance, i.e., the length ol the minimal path that connects 
x and y. 

Theorem 4.1. Assume that the graph (A, £) is 2D-like, i.e., there exists a constant C such 
that for all x € A and all integers k, 

#{yeA:d(x,y) = k}<Ck. 

Then for all u 6 [0,1], (3 G [0, oo), 5 G |N, there exists a constant K that depends on the 
graph only through C, independent of x,y, such that 

{SlSl) = \S(S + 1)F(E^ ) < f 

6 v l °s( d ( x >y) + 1 ) 
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We prove Theorem |4.1| with the help of the following lemma, that allows to compare 
expectations with respect to two Gibbs states, and something that is almost the relative 
entropy. 

Lemma 4.2. Let A,H,H' be hermitian matrices such that Tr er H = Tr cr H = 1. Then 
for any s > 0, we have 



1 



2 n s\\A\\ 



TrAe~ H -TrAe' 11 < -Tr (H 1 - H) G - R + sA + S pi| 2 e 

s 

Proof. We start with the Taylor series with remainder; there exists r](s) £ [0, s] such that 

Trc -H+sA = Tre -H +sTyAc -H , 1„2 

where F(r)) is the Duhamel two-point function 



F(r)) = [ TrAe-*"-^ Ae-^-W"-^ dt. 
Jo 



(4.1) 



Then 



TrAe- H - Tr Ac~ H ' =-Tr( C - H+sA - c - H '+ sA ) - *- {F{ V (s) - F'(r/(s))) . 



(4.2) 



(4.3) 



The remainder can be estimated using the convexity of the integrand in ( |4.2[ ) (as function 
of t) and the minimax principle: 

< F(rj) < TrA 2 e- H+r > A < \\A\\ 2 Tv e - H +^ A < \\A\\ 2 c'""- 4 " . 



Applying Klein inequality to the middle term of (4.3 1, we get the lemma. 



Proof of Theorem \4-l\ We work in the quantum setting. Let <f> z be angles such that 
and 4> y = 0. We consider the unitary 



u = n j*'**' 



z£A 



Since U*S%U=-S% and U*S*U = we have 



(S x S y ) H ( U ) — (S X S y ) ut , H ( U ) 



Then 



We have 

= e -W.-^)# ( uTz2 , + (i _ u)Qzz ,) 

= uT zz , + (1 - u)Q zz , + (cb z - c/> z ,)[S 2 z ,,uT zz > + (1 - u)Q zz ,] + 0{{<j> z - <j) z ,) 2 ). 



(4.4) 
□ 
= 7T 

(4.5) 

(4.6) 
(4.7) 

(4.8) 



The first equality follows from rotation invariance, see Lemma 3.1 Combining (4.7) with 
Lemma 4.2 and the equation above, we get, for all s £ (0, 1], 
1 



(s x s y ) < sZb 



,. Y £ tt z -^)Tr[S 2 z ,,uT zz , + (l-u)Q zz ,]e-e< U) +^ 



Ci 
s 



(4.9) 



{z,z'}e£ 

Here, C±, C 2 are constants that depend on (3 and S, but they do not depend on the graph, 
nor on s £ (0, 1] and x, y £ A. The term with G\ comes from a double commutator, 



RANDOM LOOP REPRESENTATIONS FOR QUANTUM SPIN SYSTEMS 



15 



and from higher order terms with multiple commutators. The first term of the right side 
vanishes because of symmetry (rotation around S 3 by angle ir). Had this symmetry not 
been available, a way out is to add another term with opposite angles, as in [371 149) . 
The following choice of <f> z is rather optimal: 

, = S i 1 - \:t\Z:l]t> ad{x,,)<d{z,y), 



I otherwise. 
It follows that \d>z — <t>z'\ < i — /jF 3 \ 1 1 \ j i \ , j , so that 

d(x,y) ^ 

MeS^ _ 02 ' )2 " ° 3 £ Ck (k\o g {d{ X ,y) + l)Y ~ log(d(x, 4 ») + l- 
We eventually obtain, for any s S (0, 1], 

?3c3\ s ^5 1 



(4.10) 



(4.11) 



< (StSi) < — — + SG 2 . (4.12) 

y s log(d(x, y) + 1 



We get the claim by choosing s = 1/ -y/log(d(a;, y) + 1). □ 

5. Occurrence of macroscopic loops in dimension d > 3 

The occurrence of a phase transition at low temperature accompanied by spontaneous 
magnetization and symmetry breaking was first established by Frohlich, Simon, and Spencer 
for the classical Heisenberg model in dimension d > 3 They introduced the method 

of infrared bounds and reflection positivity. The difficult extension to quantum systems 
was done by Dyson, Lieb, and Simon for the Heisenberg model with antiferromagnetic 
interactions and large enough dimension (or spin) [19] . Following observations of Neves and 
Perez [50], the result was improved by Kennedy, Lieb, and Shastry to all S £ |N and all 
d > 3 [H]. Notice that these results cannot hold in d < 2 as they would be contradicting 
Theorem |4.1| (Mermin- Wagner) . (Symmetry breaking can take place in the ground state, 
i.e., in the limit of zero temperature, see |42j and references therein.) All these results are 
proved using the method of reflection positivity and infrared bounds that was developed in 
[26] [P9] [22] [23] . See 0] for recent advances. We recommend the Prague notes of Toth 
and Biskup for excellent introductions to the subject [N21 [IT] , see also [N3j . The Vienna 
lectures of Frohlich offer an impressive glimpse of the background and of the context |N1] . 

In this section we apply the method of infrared bounds and reflection positivity to the 
model of random loops. We introduce a model with external fields and cast its partition 
function in a reflection positive form. We do not use the methods of [TH |4T] directly, but 
these articles are lighting the way. 

5.1. Setting and results. Let denote the cubic box in Z d with side length L and 
periodic boundary conditions, and let (Al,£l) the graph where edges are nearest- neighbors. 
If a formal definition is needed, consider the quotient set = (1*/LZ) d . The Euclidean 
distance in 1 d has a natural extension in this set, and we define £l as the set of unordered 
pairs {x,y} C A^ such that the distance between x and y is 1. In the sequel, we identify 
Al with the set 

A L = {x e 1 d : -f < Xi < § , i = 1, . . . , d} (5.1) 
Let us introduce a notation for the loop correlation between (0,0) and (cc, i) : 

K,(x,t) =¥(E 0>X)t ). (5.2) 
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The Fourier transform plays an vital role and we use the following conventions: 



«(fc>T)= E / 



13 (5-3) 

dte- ikx - iTt K(x,t), 



xEA L 

where k belongs to the set 

A* = {k e : -7T < hi < ir,i = 1, . . . , d}, (5.4) 
and re ^-Z d . It is useful to write down the inverse transforms: 



fceA* 



(5.5) 



/9L a 

Let e(/c) denote the "dispersion relation" of the discrete Laplacian, 

d 

e(k) = 2\J(1 — cosfcj). (5.6) 
i=i 

Let us note that the Fourier transform of k is related to the expectation of the lengths of 
macroscopic loops: 



K @) = ^S{^) = ^o ( o). 

In order to state the main theorem, we need to introduce certain integrals: 



(5.7) 



(5.8) 



An) 1 f / , n ^e(fc + 7r) 

= ,. / ( u+ 1- u) — ( 



Here, e(fc + 7r) = 2 + cosfcj), and (•) + denotes the positive part. 
Theorem 5.1. Let d > 3 and u € [0, |]. We have the two lower bounds 



"i- g ^4"V^o,e 1 , o); 

The proof of T heore m |5.1| can be found at the end of Section |5.3| 
Using Theorem 



1 fl_ 2S±1 

lim lim — = V P Bn „ > I ^ 



3.4 



(b) we obtain a lower bound for E ( -j^r ) , so that a positive lower 
bound implies the occurrence of macroscopic loops at low enough temperatures. The claim 
also holds for d = 2, but with the order of the limits over f3 and L interchanged. This 
theorem is proved in Section 5.3 The lower bounds involve P(i?o,ei,o)> which, curiously 
enough, needs to be small in the first bound and large in the second bound. We therefore 
get a positive lower bound whenever 

y / P(iW) > 2 -^I { d u) or ^(^o)< - (5. 9) 
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At least one of the two inequalities above holds when 



25+1 t-(w) , 
v^2 J d < 



(2S+1)J^ 



(«) r(tt)\-i 1 _ q(u) 

J d ) 2= b d ■ 



(5.10) 



Values of /i and J d " ; can be found numerically. The most relevant values are listed on 
Table [lj Occurrence of macroscopic loops (at temperatures low enough) follow for S = \ 
and d > 3, and S = 1 and d> 5. We shall estimate P(i?o,ei.o) from below in Section [7] in 
the case u = 0, see Lemma |7.4[ which implies macroscopic loops in d = 4 as well. They 
should occur in d = 3 as well. 



r(«) 



d 


r(0) 
J d 


t(0) 


o(0) 
^d 


r(|) 
1 d 


T (i) 

J d 


o(|) 

s d 


2 


0.646803 


1.3932 


0.24487 


0.488639 


1.28576 


0.39209 


3 


0.349882 


1.15672 


0.61150 


0.278359 


1.11536 


0.76904 


4 


0.25395 


1.09441 


0.84128 


0.208291 


1.06995 


0.99786 


5 


0.206875 


1.06754 


1.00465 


0.173274 


1.05021 


1.15760 


6 


0.177712 


1.05274 


1.13480 


0.151065 


1.03928 


1.28458 



Table 1. 
and u = 0, 



Numerical values of the integrals of Eqs (5. 
l 

2 ' 



for d 



..,6 



5.2. Reflection positivity of the model of random loops. The setting needs to be 
modified somewhat in order to reach a reflection positive form. Recall that a — (a x t) 
denotes a space-time spin configuration where the possible values are the eigenvalues of spin 
operators: a xt E {—S, —S + 1, . . . , S}. We will need a more symmetric index set for S ^ \, 
and we therefore consider the regular simplex T2S+1 that contains 2S + 1 elements. We 
embed 72S+1 in WL 2S in such a way that the points are given by vectors a, a £ {— S, . . . , 5*}, 
that satisfy 



||a-6|| = 1 if a f b; 
_ 1 



2(25+ r 
S 



m\ 

s 



if a b; 

(5.11) 



2S + 1' 

a • b — for any fixed a g {— S, . . . , S}. 



b=-S 

These properties are straightforward if we consider the projection of the unit vectors of 
jj2S+i on £ ^.jjg hyperplane perpendicular to (1,1,..., 1). More precisely, we define 

d= ±(e a - 2^i(l,l, ■■■,!)), (5.12) 

where e a is the unit vector in M 2S+1 that corresponds to a. Since all vectors a belong to a 
common 25-dimensional subspace, we can view them as elements of R' 2S . 

We denote ay' the ith component of the vector a. Let v be a real "field", i.e., a function 
— > R. The key object is the following partition function Z(v): 

z(v) = I d PL (0 ]T cx p{ E / d *[^S " s v% v - - v y ] ~ K u * - v v) 2 



(5.13) 
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Notice that Z(0) = Z^(/3,A). 

Let i = 1, . . . , d and £ = |, |, . . . , L — ~. We let R^g denote the reflection A —> A across 
the edges {x, y} e £ with Xi = £ — \ 1 y<i = £ + \ . That is, 

Ri,e{xx, ... ,xu ... ,x d ) = {x\, ... ,-(xi -£),... ,x<i), (s.w) 

where — (xi — £) is taken modulo L. Let A' 1 ) be the set of sites "to the left" of the reflexive 
plane, and A^ the set of sites "to its right" . Namely, 

AW ={ieA:i,=f-|L,..,f-|}, 

A (2) ={xe A:x t = £+\,,..,l+\L}. 

We write v = (v^',v^) where «W = is the restriction of v on A™, and i?i>W is the 

field on A' 2 ) such that 

(Rv^) x = (5.16) 

for any a; € AW. same for u( 2 '. We can now formulate the property of reflection positivity. 

Proposition 5.2. Assume u £ [0, |] and = 2, 3, 4, . . . for any reflection R, we have 

Z( W «i, (2) ) 2 < Z(v^,Rv^)Z(Rv^,v^). 

Proof. Let v x = v x e\, where e\ is the first unit vector in R . We rearrange the partition 
function (5.131 as follows. 



Z(v) = J ' dp L (£) ^ ex p|~ E J dt (° xt + \v x - o yt - \vyf 



<TSS(C) k {x,y}££' 



(5.17) 



{x,y}e£~ 

In order to reach a reflection positive form we split the Poisson process. Let the intensity l' 
be defined by 



(5.18) 



(5.19) 



, ( ( b a a b b b a a \ \ 

= }) =l-2u if a ^b, 
and = for all other subsets A of {—S, . . . , S} 4 . Let t" be defined by 

V I a a J ae{-s : ...,s}/ 
and = on all other subsets A. Using Lemma 2.2 one verifies that 

z( V )= fdpAt') /dpxn X! cxp {~ E f toi?*t + ¥ 

Next, observe that given a realization £' of p t /, we have 

/"dp tW (n F ^)= E F(„)[dp L „(n n ^ 



(5.20) 
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This holds because space-time spin configurations are constant at the transitions of 
Consequently, given £' and a £ £(£'), 



xx xxy 1 N + N °^,t,<Ty.t 
{z,y}e£t=l 



fdpAn n n n( 

(a:,I/,t)e«" {x l3 ,}e£:t=l 

=a)Too n n(i-^(i-^-., )) 

{z,y}e£t=l 

= expj- ^ ^ dt(l-^ xt ^ Bt )|. 



{z,y}e£' 



(5.22) 



We now use 1 — 5 a U<T t — (a x _ t — &y,t) 2 and we obtain 



Z(v) = f dp L >(g) ^2 ex p{- X! / &t{&xt + \v x -a 



_ 1 ^ \ 2 



(5.23) 



The measure p t / is nonnegative if u <E [0, |], and it is reflection-symmetric. See Fig. ^ for an 
illustration. 



02Q 



12Q 



Q02 



6 12 



Figure 4. Reflection in the space direction. We chose S = 2 in this illus- 
tration. The white circles mean that the spin must flip from one prescribed 
value to the other, independently on both sides. It follows that the Poisson 
point process p L > is reflection symmetric since the constraints are identical 
in both halves. 



Let us partition £ — £( 1 ' U £^ U £ , where £W and £^ are the sets of edges of and 
A^ 2 ), respectively, and £ is the set of edges with endpoints in both A' 1 ) and A^ 2 ^. We can 
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write 
Z(v)- 



£x[0,/3] 



dp AO 



/ <W(£ (1) ) V exp{- V / dt(a xt + \v x - a yt - \v v ) 2 \ 

f 13 2} 

CX P\ X! J dt(a xt + \v x - a yt - \v y ) j. 



1 <-> 2 



(5.24) 

Here, £ A (i)(£( 1 - ) ,<£) is the set of space-time configurations on x [0,/3] that are compatible 
with the specifications of £W and £. 

The couplings between both sides can be handled with the help of extra fields, as in [26] . 
Namely, we have 

expj- ^2 J dt(a xt + \v x - a yt - ^v y ) 2 ^ 



{x,y}<££ 



N 



j™ n n ex p{- n fa* + - s vt - 1%) } 

N^s 



(5.25) 



Hm n n 



=i 



do? 



x,y,t 



(— V 

V 4tt/3/ 



Let us introduce 



exp 



{x,y}e£W J ° ' L {x,y}e£t=l 



C*x,y,t ■ {<?x,t + \vx) 



(5.26) 



With this definition, the partition function Z(v) takes the form 

JV 



Z(v) = lim 

N-Hx 



Fn(v {1 \£, {ctx, y ,t})F N (v( 2 \£, {a x ,j/,t})- (5.27) 

Using the Cauchy-Schwarz inequality, and retracing our steps backwards, we obtain the 
claim of the proposition. □ 

Proposition 5.3. Assume that u 6 [0, |] and 9 — 2,3,... Then Z(v) is maximized by 
v = 0. 

Proof. Existence of maximizers is clear, since Z(v) is continuous and positive, and it tends 
to as any v x is sent to infinity. Beside, a constant field is equivalent to the zero field. 
The argument is then standard, one uses reflections to construct a sequence of maximizers 
where more and more fields are constant, until they are all identical. See Fig. [5] for a quick 
illustration, and the references Q15] or the notes |N21 IN3| for more details. □ 
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(HHH» 



m 



Figure 5. Starting with a maximizer, reflections yield further maximizers 
where more and more values are identical. 



5.3. Infrared bound for the correlation function. The name "infrared bound" refers 
to estimates of the Fourier transform of P(-Eo,x,o) around k = 0, that is, for low "frequen- 
cies" . We first use Corollary |5.3| in order to derive a bound on the Fourier transform of 
the Duhamel function. Then we use Falk-Bruch inequality to extend it to the ordinary 
correlation function. These steps follow [19] ■ 
Recall the definition of k in Eq. (5.3). 

Proposition 5.4. Assume that u £ [0, h] and 9 = 2, 3, . . . For k E A* \ {0}, we have 

**,»>< 

From now on, we let (v,v') denote the inner product in £ 2 (A), that is, 



Proof. We have 



E(o- ,o • o-x,tj- 



It follows that 



S 

S 



2S . - 

x£A 

Let A denote the discrete Laplacian, such that 



r«(fc,0) = V / dfe^Efoo-SW). 
1 ...Jo 



y.{x,y}eS 



We have 



(5.28) 



(5.29) 



(5.30) 



(5.31) 



(5.32) 
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We choose v — cos(fccc). Expanding around v — to second order, using — Av = e(k)v, we 
get 



0(ry 4 ) 



1 + i77 2 /3e(fc) 2 / di e((?$,v)(c?<J,i;)) - Wps{k) 



Z(0) 



We calculate the expectation: 

EUa^Q ,v)(a^t ,vfj = ^ cos fcrcos ky Kuf^a^ 



(5.33) 



= cos /cx cos /c(x — z) E^a^a 

I,2£A 



(1)^(1) 

z.t 



(5.34) 



The last line is obtained by replacing y by —y, by using translation invariance, and with the 
substitution z = x + y. Observe now that 



We now get the Fourier transform of correlations: 



k(0, z, t). 



(5.35) 



zGA 



cos fcx K,(k, t). 



The last identity holds because n[k, t) is real due to lattice symmetries. Then 



1 


2(25 4 


-1) 


1 




2(25 4 


-1) 


1 




2(25 4 


-1) 



z£A 



(5.36) 



Z( V v) = Z(0)(v,v) 



1 + 4PS+I) V 2 Hk) 2 / dt /£(*, t) - \r?fe{k) 



The bracket is negative for small rj, so that 



1 



25 4-1 



e{k) / dt1i(k,t) < 1. 



The proposition follows. 



(5.37) 

(5. 38) 
□ 



The next step is to transfer this bound to the Fourier transform /?(fc,0), see Eq. (5.3). 
We need to compute the double commutator of the Falk-Bruch inequality. Let us introduce 
the operators Q^ y and Tj^ in H x , v by 



(a, b\Ql 6 y \c, d) = {a- c) 2 5 ab 5 cd , 
(a,b\T™\c,d) = (a-c) 2 S ad S bc . 



(5.39) 



The peculiar notation is motivated by the fact that these operators are given by double 
commutators with Sr. 
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Lemma 5.5. 

Qxy = ~[^u [Qxy, S x ]] — [S'y , [Qxy> S x ]], 

7^33 r r*3 \rp c3]i r c*3 \rp q3]] / c3 c3\2ti 

1 xy — [ D xi l-Lxy, ^> x \\ — l°y V-xyi °xJJ — K°x °y) 1 xy 

The proof of the lemma is immediate by looking at the action of these operators on basis 
elements. Next we consider 

£3 = J2e- ik *Sl (5.40) 

and we introduce 

fit) _ /rp33 



(5.41) 



(u) 

where the expectations arc taken in the Gibbs state with Hamiltonian . These number 



(5.42) 



turn out to be related to the probabilities of Eq , see Lemma 5.8 below. 
Lemma 5.6. 

([S 3 _ k , [H^,S 3 k }\) = |A| (ue(k)r[ u) + (1 - u)e(k + tt)t (u) ). 

Proof. We perform the computations for and separately, and we use linearity to 
get the result for ' , 

[H£\si]=Y,*- ikx K^s 

= - E e- ik *[T xy ,Sl) 

x,y-{x,y}e£ 

= J2 e- ikx T xy (S 3 x -S 3 y ). 
x,y-{x,y}ee 

[S 3 _ k , [flf ,Sl}} = - £ e~ ik * [ e ifc * S 3 X + S 3 , T xy (S 3 X - S 3 )} 
x,y-{x,y}e£ 

— [j x -r e o y , i xy (D x D y )\ 

x,y-{x,y}e£ 

= J2 (l-e- iHx - y) )(S 3 x -S 3 fT xy 
= 2 E (l-cosfcCx-y))^-^) 2 ^. 

{x,y}€£ 

Next, we perform the calculations for . For the commutator, we get 

{Hl°\sl} = Y,e- ikx l<\s 3 x } = - E e- ikx [Q xy ,S 3 x ]. (5.44, 

And for the double commutator, we get 

[S 3 _ k ,[H<£\sl}}=- J2 [S 3 x+e- ik ^S 3 ,Q xy (S 3 ^S 3 )}. (5.45) 

x,y{x,y}e£ 

The result follows from Lemma IB~5l □ 



(5.43) 
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Proposition 5.7. 

WlKk) < ^ S(S+1 l (2S+1) ^ + (1 - u)r^ + S{S + 1} . 

Proof. We use the inequality (3.261 with 

(A*A + AA*) = (S 3 _ k Sl + SlS 3 _ k ) = 2\A\(S$s!)(k). 

We have 

(5|,5|) Duh = |S(S + 1)|A|k(M). 

Using Lemma [5T6| we get 



(5.46) 



(5.47) 



2|A|(sfi|)(fc) < ^5(5 + 1) | A|«(fe, 0) 7|A| (^(fc)^ + (1 - + tt)^) 

25(5+1) 



+ 



3/? 



The claim now follows from Proposition |5. 4 



|A|7c(fc, 0). (5.48) 

□ 



The next lemma allows to get upper bounds for Tq 11 ^ and t± . 
Lemma 5.8. For any u G [0, 1], we have 

(a) T^ = IS(S + l)(2S + l)F(E^ eu0 ). 

(b) T^§5(5 + l)(25 + l)P(£ + e > 

Proof. We start with (a). Given a realization uj of the Poisson point process p Ul we denote 
u) the realization where we add a double bar on the edge (0, ei) at time t — 0. We have 



Z(")(/3,A) 



/ dp„(w) £ (o- ,o+ - cr , -) 2 

•*'- B §,ei,0 ff££(ffi) 



(5.49) 



The first integral is over realizations u; such that (0,0) and (ei,0) belong to different loops. 
The additional double bar merges these two loops, so that ui G £b,ei,o and cr 0! o+ = o"o,o- f° r 
all compatible configurations. The first integral is then zero. 

The second integral is over u> G E 0ei0 . In this case, Q G E 0ei o if oj G and 

w G Eq Ci if w G Eq Bi0 . We get in the first case, and 

£ Ko + -^o-) 2 = (2^ + l)' £ ^'- 1 £ (a-fo) 2 

a-eE(cfl)) o,6=-S 

= (25 + §5(5 + 1)(25 + 1) 

in the second case. We used the identity 

^ ]T (a-6) 2 = |5(5+l)(25+l) 



(5.50) 



25 



(5.51) 



a,b=-S 



which holds for all S G |N. 
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The proof of (b) is very similar. We have 



» 



Z(")(/3,A) 



dp u (u)l Eoi * li0 (w)lEl eii0 (&) E ( CT 0,0+ -°0fi-) 2 - (5-52) 



It differs from (5.49) in that Q contains an extra cross instead of a double bar. The two 
indicators are equal to 1 E + and we get the claim of the lemma. L 



0,ii-i ,0 



We can now prove the claim about the occurrence of macroscopic loops. 



Proof of Theorem 5. 1 We have 



((So) 2 ) = zM \^ K) J dp«M E CT oo = \S{S + 1) 



Then 



hS(S + 1) — (S S X ) 



x=o |A| 

Using the infrared bound of Proposition |5.7| we obtain 



1 (s§sfxo) + 1 l E OOTX*)- 

1 1 fceA*\{o} 



(5.53) 



(5.54) 



i , /gfg+iKgg+i) 



|A| 



|A| 



fceA*\{0} V 
S(S+1)(2S+1) 



( M ) g(fc + 7r) 

e(k) 



1 \- 1 

/3IAI ^ eikY (5 ' 55) 



The first lower bound of Theorem 5.1 follows immediately from Eq. (5.551, Lemma 5.8 



and the definition (5.8) of . The second lower bound is obtained as in jJT], using 
\- E (S$SlHk)cosh (5f^)(fc)(e ifcl + e" ifcl ) 



|A| 



2 I A I „ , 
h({S 3 Sl) + {SlS* ei }) 

§S(S+l)P(£ , ei ,o). 



(5.56) 



It follows that 



5(5 + l)P( J BQ, ei ,o) = ^T^)(0) + ^| E (^|>(fc)E C0S ^ 



fceA*\{o} 



< ^ ( 535i) ( o, + ij^¥^^ j; ^« +(1 -„ )T «£a±E)(|:cos t( ) + 



S(S+1)(2S+1) 



— y — ■ 



(5.57) 

The second lower bound of Theorem 15.11 now follows from Lemma 15.81 and the definition 
Qof^"'. □ 
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6. Reflection positivity in space and time 

This section describes an extension of the method of reflection positivity and infrared 
bounds to the space-time. As it turns out, the results are slightly less sharp than in the 
former section. It seems nonetheless useful to include this for several reasons. First, the 
idea is natural and it may save the readers some efforts. Second, it might prove useful in the 
future. It would indeed provide better bounds if reflection positivity could be extended to 
U > | . Reflection positivity in space and time has been independently (and indeed, earlier) 
proposed by Bjornberg for the quantum Ising model in order to prove interesting results 
about critical exponents in d > 3 [13]. The method described here shares many similarities. 

The main result is the following. It is very similar to Theorem |5.1| it is actually identical 
when u — \, but it is not as good when a£ [0, h). 



Theorem 6.1. Let d>Z and u G [0, §]. We have the two lower bounds 
lim lim E(^) > 



£(o,o) n > jl-(25+l)J^V(l- M )P(£ , ei , ); 
P^Ll^V ) - \p(E , ei , ) - (2S + l)I i d i) y/(l-u)¥(E 0iei ,o). 



The numbers 1^ and are given by the integrals in (5.8) with u = ~ and their values 
for d = 2, . . . , 6 can be found on Table [T] 

The rest of the section is devoted to proving this theorem. The method of proof does 
not rely on analytic inequalities such as Falk-Bruch, and it will be more attractive to some 
readers. Its difficulty is about the same, though. 

We generalize the notion of partition function with external fields, by considering fields 
v : A x [0, /?] per - > K- That is, v also depends on the "time" parameter. We prove that 
the partition function has a local maximum at v = 0, and we obtain a generalized infrared 
bound for the Fourier transform in space and time of the correlation function. 

Let V Co be the set of fields v : A x [0,/3] per — > K where v xt is twice differentiable with 
respect to t, and | | < Co for every x,t. We introduce the partition function by 



Z(v) = [ dp u (u>) ^ CX p{ [ dt fix} ~ &y})( V x,t ~ Vy,t) ~ \(v x ,t ~ V y ,tf 

rr ez A u ' 



(6.1) 



The constants a and b will be chosen later. 

Proposition 6.2. Assume that u E [0, i]. For every v E V Co , there exists v* = (uj) € V Co 
that depends on t but not on x, such that Z(v) < Z(v*). 

The proof can be done by extending Proposition |5.2| to the partition function above, and 



by repeating the proof of Proposition 5.3 It turns out that the time-dependence of the fields 



and the extra terms do not play any role. 

Proposition 6.3. Assume that u E [0, |] and that b > 2rfa 2 (l — u)«;(ei,0). Then there 
exists cq > such that Z(v) < Z(0) for every v E V Co . 



Proof. Because of Proposition |6.2| it is enough to prove it for space-invariant fields only. We 
must show that maximizers are time-invariant as well. We have 

Z(v) = lim Zpf(v) (6.2) 

JV-S.OO 
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where 
z n{v) = / dp u (u) 



E 



E 

te#{i,...,iV} 



a(<? (1) „ 



CTw)K + 4-W t )-^(w t+ A-Wt) 



(6.3) 

Indeed, we have discretized ^Sr and 4^ and used the discrete integration by parts. We 
take the limit along N e 2N. 

We apply "horizontal reflections" across the planes determined by t — rij?, n = 1, . . . , N . 
Then Z N (v^, ) < Z N (v {1 \ Rv^)Z N {Rv^\ w (2) ) from the Cauchy-Schwarz inequality 
— the situation is actually simpler than in the proof of Proposition |5.2| as we do not need to 
introduce extra fields to decouple the two parts. See Fig. [6] for an illustration. Proceeding 
as before, we obtain the existence of a maximizer v* for Zn with jigsaw shape, namely 



for some constant |c| < Co/3/2. 



, ,« c 

( — 1) ~3~ — 



(6.4) 



0/N 



Figure 6. Reflection in the time direction, across the dotted lines. There 
are two dotted lines because of periodicity. 



We need to show that c = 0. We have 

46c 2 1 A | 2ac 



Zn(v*) = / dp u (uj) ^2 ex P] 



crGS(w) 



(6.5) 



We need to show that the term linear in c is actually quadratic. This can be done using 
the following trick. When integrating over realizations lo with the measure p Ul we replace 
each transition (x, y, t) G cj by h transition at t, and | transition at t + j?. Because of the 
alternating sign in the maximizer v* , we obtain 



Z N (v*) = e—F 



-|A| 



/^.HE n hU{ 2 -f(t ( x 

ff£u (x,y,t)e.u \ ^ 



exp 



2ac , (1) 
' 1 *.t+# 



? (D 



.(1) 




- ?s 


*s - 


*<>}) 





x,i u y,t, 



(6.6) 
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The correction 0(]y) is due to the realizations w where transitions occur at almost the same 
location and time. 



If the transition is a cross we have er 



a Vt t and a y t+ p_ = a Xt t and the corresponding 



factor is 1. If the transition is a double bar, we get the factor 

_;8a 2 c 2 

/3 v "=M+# 



cosh 



,i r 8a 2 c 2 f _(i) ^(1)^2 . 4 . \ 

)J = ex P|^2-(<; +# -<0 +0(c )}. 



(6.7) 



Let T(w) denote the set of double bars that are present in the realization lu. We then obtain 
Z N (v*)=e-^ |d Pu H^exp{^ J2 ( (? S + A"^ ) ) 2 + 0( C 4 )}+0(1). 



(i,v,{)eT(u) 



(6.8) 

Notice that 0(c 4 ) is bounded uniformly in N, and O(^) is bounded uniformly in c. They 
depend on |A| and P, on the other hand, but it does not matter. Let A C O be the event 
where a double bar occurs on the edge {0, e{\ in the time interval [0, -ft]. Expanding the 
exponential, and using translation invariance in space and time, we get 

Z N ( V *) = Z N (0)[1-^\A\]+^\£\N [ d Pu (u) J2 «i-^) 2 + 0(cVO(|). 

(6.9) 

We now estimate the contribution of the realizations with an enforced double bar on (0, e\) 
at a time close to 0. As before, we denote u> the realization lj with an extra double bar at 
{0, ei }x0. 

cr££(cj) 



lim , 

N^oo Zjv(0) P J A 



(1 — u) lim — 



(6.10) 



= (l-u)P(£ - eii0 ) 
< (l-«)P(^ ieii0 ). 

The sum over space-time spin configurations that are compatible with uj G Bo.ei.Oj wa s 
computed as follows: 

E (^i--S) 2 = (25 + l) l ^ )| - 1 £ (a«-6W) 2 



a.b=-S 
2S 



(Mf + 1) |/*-)M^£ ^ (3»-SW) a 



(25 + 1 



|£( w )l-iJ_ 

25 ^ 



»=l ,6=-S 
S 



(6.11) 



a,b=-S 



= (25 + l) |£(w)l . 



We used Eqs (|5.11|). We have obtained 

1 



Z N {v*) < Z N (0) 



Abe* 



|A| + ^d|f|(l-«)/s(ei,0) +0(c i ) + 0(j r ). (6.i2) 



We see that c = is local maximizer whenever ^ > ^§-^(1 — u)fc(ei,0), which yields the 



relation between a and b that is stated in the proposition. 



□ 
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We now use Proposition 6.3 in order to get a generalized infrared bound for n(k, t). This 
is similar to Proposition |5.4| 



Proposition 6.4. Assume that there exists cq > such that the partition function in Eq. 

or every v € V Co . T 

K(k,r) < (25+1)- 



(6.1) satisfies Z(v) < Z(0) for every v € V Co . Then for all (fc, r) ^ (0,0), we have 

e(fc) + 46r 2 



(e{k) + aT 2 y 



Proof. We choose v Xyt — cos(fcx + rt). We have 

rP 



Z(qv) = J dp u (uj) ^ ex P 



dt 



T)$$,Av. lt ) - lr] 2 (v. tt ,Av., t ) 



E / 



^(l)C «a:,t , 2 



dt 2 



.13) 



We now use — Aw = g:(fc)w and -pf = t 2 u, and we get 

Z(nv) = [ d Pu (u) ex p( f P dt\r,(e(k) + ar 2 ) (a [ \\v, t ) - +46r a ) («.,*,«.,*) 

r/3 r /3 



,t' ' t/ -,t' 



y^w) E (l + h 2 (eW+ar 2 ) 2 J dt dt'(a^,v. tt ) (c? ( 

- \<rf{e{k) + 4br 2 ) J dt(v. tt , v., t ) + 0(r/ 4 )) . 



cr££(w) 



We have 



(ct. j , u. :4 ) (07^, w. ; t') = ^ cos(kx + rt) cos(ky + Tt')a x ' t a, 

x,y£A 



^^(1)^(1) 
V,*'' 



(6.14) 



(6.15) 



and, using the symmetries of the cubic box and Eqs (5.11 1, 
1 

Z(0) 



(7€S(u) 



z (°) K-..^ 
1 

2(2S*+ 1) 



dp„H(25 + l)l £ ^l- 1 ^ £ ||a|| 
K(y — x,t' — t). 



a=-S 



(6.16) 



Finally, we obtain the Fourier transform of correlation functions: 

rP 



rP rP 

/ dt dt' cos(fcx + rt) cos(ky + Tt')n(y — x, t! — t) 
Jo Jo x,yeA 

rP rP 

= dt dt" ^2 c °s(fca; + rt) cos(k(x + z) + r(t + t")n{z, t") 

rP rP 
= / V cos(fca; + ri)Re e ikx+iTt / dt" e ikz+wt n(z,t") 
Jo xeA Jo 
rP 

= / dt >J cos(fca; + rt) 2 n(~k, — r). 



(6.17) 
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We actually have «(— k, — r) = n{k,r) because of symmetries. We have obtained 



Z(r)v) = Z(0) 



4(25- 



— V(e(fc) + ar 2 ) n{k,r) / / di(w. )t) r;., t ) 



r; 2 (e(fc) + 46r 2 ) 2 / di(v., t , u.,t) + 0(ry 4 ) 



I„2 

4 



The conclusion follows. 



(6.18) 



□ 



Corollary 6.5. Assume u £ [0, |]. Then for all (fc,r) 7^ (0, 0) 7 we /icwe 



/e(fe,r) < 



25 + 1 



e(fc) + 



8d(l-«)«(ei,0) 



Proof. It follows from Proposition |6.3| that the bound of Proposition |6.4| holds with b = 
2da 2 {\ — u)n{e\,Q), for any a > 0. We get the result by optimizing over a. □ 



Corollary 6.6. Assume u €E [0, |]. Then for all k^0, we ho 



8(*,0) = i E «(*, r) < (25 + 1} " )K(ei ' 0) coth(/3v/ 2rf(l - ^R^oRfc)) . 



This corollary follows from the identity X^nez w 2 +g 2 = f c °th(7r£). We can now prove 
Theorem 16.11 



Proof of Theorem 6.1 For the first bound, we use 



1 = k(o, 0) = E E «M 



|A|/3 



re3fZ\{0} 



|A|/3 



fceA*\{o> re- ; 



(6.19) 



The first term is equal to E( ) . The middle term vanishes in the limit |A| — > 00. The last 



term can be bounded by Corollary |6.6| recalling that k is real because of lattice symmetries. 
We get 



lim E, 

|A|-hx> V f3\A\ 



coth(/3 v /2d(l - u)K(e ll 0)e(k)) 



(27T)° 



dk. (6.20) 



The cotangent disappears in the li mit f3 — > 00 by dominated convergence (in d > 3) and we 
obtain the first bound of Theorem 6.1 



since the integral of l/y/e(k) is equal to J^ 2 ' j\j2d. 
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For the second bound we use the sum rule of [STj . Using invariance under lattice rotations, 
we have 

1 _ d 



1 IM 1 IP re^z\{o} 1 IP fceA*\{o}re^z; i=i 



(6.21) 

As before, the first term is equal to E( p?£? ) and the middle term vanishes in the limit 



|A| — > oo. Using Corollary 6.6 we get 
- L 

lAjToo^V /3|A| 



lim E (-Jrr) ^ K ( e li °) - ( 25 + !) V 72 ^ 1 - w)«(ei, 0) 



x * / coth(V2d(l-u)« (ei ,0) £ (fc))^ x dfc ^ 

d(2TTY y^jfe) 7 + 

The cotangent again disappears in the limit (3 — > oo and the integral is equal to I\ 2 /y2d. □ 

7. Specific models of interest 

It is time to give flesh to the quantum models under study. The Hamiltonians were 
defined in terms of operators T xy , Q xy , and P xy , that were chosen for their mathematical 
convenience rather than their physical relevance. In this section we discuss the special cases 
S = ^ and S = 1 in details. 

7.1. Spin-i models. Let us start with the spin-^ Heisenberg ferromagnet, the model 
that was studied by Conlon and Solovej |16j , and by Toth who introduced the representation 
with "crosses" [S3]. The parameters are S = \ and u = 1. The loop parameter is thus 9 = 2. 
Using S x ■ S y — \T xy — j Id, we find that 

E (S X -Sy+\). (7.1) 

{x,y}££ 

Since only transpositions are present, the loop representation can be seen as a model of 
"spatial permutations", i.e., bijections A — > A. Each loop corresponds to a permutation 
cycle, and the length of the loop is equal to /3 times the number of vertices in the cycle. 
Next, we discuss the spin-^ Heisenberg antiferromagnet. Let S = \ and u = 0, and 

consider the Hamiltonian = — Y^{ x y} ^ x v We nave S x ■ S y = |ld — \Pxyi 80 that 

ffi 0) =2 Y, {S X -Sy-\). (7.2) 

{x,y}££ 

This is indeed the Hamiltonian of the Heisenberg antiferromagnet. We cannot use Theorem 
|3.5| because the spin is half-integer. But if we assume in addition that the graph is bipartite, 
i.e., A = Aa U Ab such that all edges of £ involve one site in Aa and one site in Ab, the 
Hamiltonian is unitarily equivalent to H^: 



4 o) =(n ^K^n 

i6A B iEA b 



e™ s * ). (7.3) 
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Then we can use the probabilistic representation. It only involves double bars since u = 0, 
and it was introduced by Aizenman and Nachtergaele 3J. Spin correlations are given by 

(s 1 x sd = \(-ir-y¥(E x , yfi ) (7 .4) 

for i = 1, 2, 3, where (— l) x ~v is equal to 1 if x, y belong to the same sublattice, —1 otherwise. 

The case of frustrated systems where the graph is not bipartite is currently attracting a 
lot of attention by condensed matter physicists. The probabilistic representation does not 
apply, because the weights of loops would carry signs. 

Another special case is the spin-| XY model, taking S — | and u—\. We can check 
that 

= ^(S^Sy — SgSy + S^y) "f" J" (7-5) 

It follows that 

Hp = -2 (SlSl+S*S* + {). (7.6) 

{x,y}E£ 

A consequence of Theorem |3.3| is that 

0< <(SlSl) = (S*Sl). (7.7, 

Neither the second inequality, nor the positivity of the latter correlations, seem to be im- 
mediate. 

An additional motivation for the XY model comes from the fact that it is equivalent 
to the hard-core Bose gas. This is well-known, see e.g. [HJ [2J, but we recall it for 
the convenience of the readers. Notice that the present XY model differs somewhat from 
standard conventions, since interactions are between spins in the directions 1 and 3 rather 
than 1 and 2, and this requires a few modifications of the usual correspondence. 

Let a x = Sl+iS^ and its adjoint a x = S x — iS x . These operators satisfy the commutation 
relations 

[a x ,a,W] = for all x ^ y, 

(7-8) 

{a x ,a x } = IdA for all ieA. 
In addition, the operator for the number of particles at x is 

n x = a x a x = S% + |ld A . (7.9) 
The Hamiltonian can be rewritten as 

iJ A 2) =- ^2 (ala y + ala x y (r.w) 

{x,y}<E£ 

The relevant correlation functions are those representing "off-diagonal long-range order" 



(ata y ) = ^Tr (S* - iS 3 x )(S y + iS 3 y ) e^ H - 

ZP (7.11) 

= 2{S 3 x S 3 y ). 

We used the identity (£*S£) = (S%S%), and the fact that {SlS]) = -(S X S 3 ) = 0. 

Existence of long-range order was proved in H] in rectangular boxes, d > 3, and low 
temperatures. Theorem |5.1| does not improve these results. But we use the random loop 
representation to give a heuristic description phase transitions and symmetry breaking in 
Section [O 
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7.2. Spin-1 rotation-invariant model. It is well-known that any two-body rotation- 
invariant interaction for S = 1 can be be written as 

H A = - ^ ^lSx ■ Sy + J 2 (S X ■ Sy^y (7.12) 

{x,y}££ 

The phase diagram of the general rotation-invariant spin-1 model is displayed in Fig. [7] The 
line J2 = corresponds to the usual Heisenberg model. Long-range order has been proved 
for the antiferromagnet when d > 3 and at low enough temperatures [THl EI] • It is possible 
to check that JiS x ■ S y is reflection positive, in the quantum sense, when Ji < and that 
J2(S X ■ S y ) 2 is reflection positive when J2 > 0. Thus H\ is definitely reflection positive in 
the quadrant J\ < 0, Ji > 2. One can obtain an infrared bound for the usual correlation 
function, which allows to extend the domain of long-range order to the domain depicted in 
Fig. [7| The domain of reflection positivity presumably extends a bit beyond, but this has 
not been shown yet. The line in the direction (—1, — ~) corresponds to the model introduced 
by Affleck, Kennedy, Lieb, and Tasaki [T] in order to study Haldane's conjecture. 
We first express T xy and P xy as linear combination of S x ■ S y and (S x ■ S y ) 2 . 




Figure 7. Phase diagram of the general spin-1 system with Hamiltonian 
— Yl{JiSx ' S y + J2{S X ■ S y ) 2 ). Neel order was proved when J x < and 
J2 = in [03 HI] ; this can be straightforwardly extended to a region with 
small J2 > 0. Antiferromagnetic behavior is expected there. The random 
loop representation allows to prove Neel order for J\ = and J2 > 
(Theorem 7.3), and to prove another form of magnetic order for < J\ < 
i J2 (Theorem 7.2). The nature of the phase transition is not clear. 
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Lemma 7.1. In the case S — 1, we have 

T X y — S X ■ Sy + (S X ■ Sy) — 1, 

Pxy (S x ' Sy*) 1. 

Proof. The claim for T xy can be verified by direct calculations of all matrix elements, using 
the identity 

^Sq; ' Sy Sy ~\~ Sj, Sy ~\~ "2S x Sy. (7.13) 

Direct calculations can also be used for P xy . However, a more elegant argument uses prop- 
erties of additions of spins (see e.g. [46]). It is well-known that the eigenvalues of (S x + S y ) 2 
are 0, 2, 6 (they are equal to J(J + 1) with J £ {0, 1,2}). This gives the eigenvalues for 
S x ■ S y : —2, —1, 1, and hence for (S x ■ S y ) 2 : 4, 1, 1. Recall that \P xy is the projector onto 
the one-dimensional eigenspace of (S x + S y ) 2 with eigenvalue 0. The identity of the lemma 
is now easily checked for any vector that belongs to the eigensubspaces. □ 



It follows from Lemma 



7.1 



that i?i can ^ e wr itten as 



(7.14) 



=- (US X ■ Sy + (S X ■ Sy) 2 ~ Id). 

{x,y}€£ 

The region of parameters where the model has the probabilistic representation (with 
positive weights of the loops) is delimited by the lines J 1 = and J 2 = J\. The ordinary 
spin-spin correlation function is given by Theorem 3.5 (a), and this leads to the following 
conjecture. 

Conjecture 1. Let («/]_, J 2 ) satisfy < J\ < J 2 . For all (3 6l and all d > 1, the correlation 
function 

(SiSl) = l[P(E+ ytt ) -E(E~ Vtt )] 
has exponential decay with respect to \\x — y\\. 

Long-range correlations are only possible if long loops are present, and the vertical ori- 
entation is quickly lost. It should be possible to prove this rigorously, although it does not 
appear to be straightforward. 

The results of Section [5] n onetheless imply a phase transition with long-range order. 
It follows from Theorem 3.5 (b) that macroscopic loops are accompanied by long-range 
correlations of (S x ) 2 . Here is the main result of this article for the model with 5 = 1. 

Theorem 7.2. Let (A,£) be the d- dimensional cubic box with periodic boundary conditions 
and even side length. Assume that < 2Ji < J 2 and that d > 5. Then there exists /3q < oo 
and c > such that 



E (((^) 2 (^) 2 )-((^) 2 )((^) 2 )) 



> C 

for all (3 > (3q, i = 1, 2, 3, a;, y G A, uniformly in the size of the system. 



This theorem is a direct consequence of Theorem 3.5 (b) and Theorem 5.1 It establishes 
the existence of a phase transition with symmetry breaking, since uniqueness of infinite- 
volume Gibbs states implies exponential decay of all correlations. Notice that similar mag- 
netic properties have been recently proved in a plaquette orbital model by Biskup and 
Kotecky [12] . 

The case u = in a bipartite graph is different. Only double bars occur in the loop repre- 
sentation, and the vertical orientation displays alternating properties. Namely, P(E x y t ) = 
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whenever x, y belong to the same sublattice, and P(E^ y t ) = whenever x, y belong to dif- 
ferent sublattices. The following result can be proved. 

Theorem 7.3. Let (A, £) be the d-dimensional cubic box with periodic boundary conditions 
and even side length. Assume that J\ = 0, J2 > 0, and that d > A. Then there exists 
0o < 00 and c > such that 

ftp E (-^-"(sisi) > c 

for all ft > fio, i = 1, 2, 3, x, y € A, uniformly in the size of the system. 

The system therefore displays Neel order at low temperatures. It is tempting to expect 
antiferromagnetic phases that are characterized by staggered magnetization. In other words, 
all pure Gibbs states should be of the form 

(")ji= lim lim < n»o n ' (7 - 15) 

U fl->0+ |A|->oo U A + h 2^x\~ 1 > u ' b " 

with 57 e § 2 . But the heuristic discussion in Section 8.2 suggests that this is not the case. 
This phase transition is not fully understood. 

For d > 5 Theorem |7.3| follows again from Theorem |5.1| We can extend it to d = 4 
by using the lower bound of the latter theorem, together with the following estimate on 
P(£o, ei ,o). 

Lemma 7.4. Under the same assumptions as in Theorem \ 7. S\ we have the lower bound 

lim P(£u ei ,o) >| + S3- 
p— ^00 



We actually prove an estimate that is valid for all S € |N, see Eq. (7.24 1 below. 

Proof. We relate this probability with the expectation of Qo. ei which can be estimated 
using a variational argument. Let (•) denote expectation with respect to the Gibbs state 
with Hamiltonian Hj^K 



1 

Z(o)"CM) 

1 



(Qo,ei) — ryinua A\ / dpoM E £o-o,o-,o-e 1 ,o- ( Wo+,o-e 1 ,o+ 

TeS(a;U(0,ei)) 

dp (w)(25 + + 1 - / dp Q (o;)(2S + 



(7.16) 



= (25 + l)P(^o, ei ,o) + ^fT P (^o c , ei ,o) 

= 2^+1 + ( 2 ^ + 1 _ ^i+l^i^^fi)- 

Using translation invariance we have (H^) — —\£\(Qo. ei ) and we obtain the identity 

W..,o) = ^(<=fp-^). 

We estimate the average energy using a trial state. Let 

D = {{x, y} e £ : x is even, y = x + ei}. (7.1s) 
We consider the dimerized state 

1 S 

|dimcr) = / 2g 1 E l a '°)' (7 - 19) 

{i,t/)£D a= — S 
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lim (H\ ) = min {ip\H^ < (dimer|ii^ |dimer). (7. 20) 



(7.21) 



Wc have 

00 

If {x, y} £ D, we have 

1 S 

(dimerlQxyldimer) = — - ^ ^ (a,a\Q x>v \b,b) =25+1 

a 1 b= — S 

If {x, y} ^ _D, we have 

1 5 

(dimer|Q X j,|dimer) = — — — - ^ (a, a, 6, fclQ^Jc, c, <i, d) = 1. (7.22) 

^ ' a,b,c,d=-S 

Since |£)| = ||A| = |, we obtain 

(dimcr| - #j 0) |dimer) = ^\S\. (7.23) 

We finally get 

r 2S+^fS + d 1 \ 

hm P(JW) > 45(5TI) l^r - 25+t)- 
The result follows with 5 = 1. □ 

8. Conclusion and outlook 

Connections between random loop models and quantum lattice systems provide many 
deep insights for each of them. The continuous symmetries of the spin systems have far- 
reaching consequences regarding the size of loops in dimensions 1 and 2. Results about long- 
range order, obtained in the quantum setting [19] [41], establish the occurrence of macroscopic 
loops when S = ^, i.e., 9 — 2. Conversely, the loop representation shows that certain 
quantum correlation functions are positive and that a phase transition takes place in the 
quantum model with 5 = 1. This phase transition had not been noticed yet. We proved it 
by adapting the method of reflection positivity and infrared bounds of [26] , rearranging the 
underlying Poisson point process of transitions so that it becomes reflection positive. 

Several of the present results should hold more generally. Long loops should be absent in 
dimension 2 for all 6 > 0, not only 9 — 2, 3, . . . Occurrence of macroscopic loops is proved 
for small S or large d (and large /3) but this could certainly be improved. The restriction 
to u < \ seems to be an inherent feature of the method; it is indeed also present in the 
quantum setting. 

Nachtergaele's extension of the loop representation for higher spins [371 SB] are different 
from the ones discussed here. Comparing the information provided by both could lead to 
new results. 

8.1. Joint distribution of the lengths of macroscopic loops. This is an intriguing 
topic for future research, both in itself and for the information it provides on the structure 
of pure Gibbs states and symmetry breaking. It seems appropriate to discuss this in details. 
Let ,L*( 2 \ . . . denote the lengths of the loops in decreasing order. Clearly, (^j, |pT, ■ • ■ ) 
is a random partition of [0, 1]. In the case of the random interchange model on the complete 
graph, i.e., 6 = 1, u = 1, and (A, £) is the complete graph, Aldous conjectured that this 
random partition has Poisson-Dirichlet(l) distribution. This was subsequently proved by 
Schramm |52) . who showed that the time evolution of the loop lengths is described by an 
effective split-merge process (or "coagulation-fragmentation"). The relevance of these ideas 
in the presence of spatial structure (that is, A C Z d ) and for 9 = 2 was explained in [25] . 
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This is backed by results for the model of "spatial random permutations" , that are rigorous 
in the annealed case [TU] and numerical in the lattice case |31j . 

We start by describing the heuristics. Our main goal is to justify Conjecture [2] below. 
We assume that the graph is regular, A C Z d , but the discussion holds more generally. 

Macroscopic loops are spread all over A and they interact between one another, and 
among themselves, in an essentially mean-field fashion. We introduce a Markov process 
for which the measure (2.3) is invariant. It is enough that it satisfies the detailed balance 
property, which can be written as follows: 

9^(udt)# c ™ in "((1 - u)dt)* douW ° bars in u R(lj,cj') 

= (udtf crosscs in »' ((1 - u)dt)* douWc bars in u 'r{w',u), (8.i) 

We have discretized the interval [0,0] with mesh dt, and R(lu,ui') denotes the rate at which 
uj' occurs when the configuration is u>. The following process satisfies the detailed balance 
property. 

• A new cross appears in the interval {x, y} x [t, dt] at rate uQ 1 / 2 dt if its appearance 
causes a loop to split; at rate u9~ 1 / 2 dt if its appearance causes two loops to merge; 
and at rate udt if its appearance does not modify the number of loops. 

• Same with double bars, but with 1 — u instead of u. 

• An existing cross and double bar is removed at rate 9 1 / 2 if its removal causes a loop 
to split; at rate 9~ 1/>2 if its removal causes two loops to merge; and at rate 1 if the 
number of loops remains constant. 

Notice that any new cross or double bar between two loops causes them to merge. When 
u = 1, any new cross within a loop causes it to split. 

Let 7,7' be two long loops of respective lengths L,L'. A new cross or double bar that 
causes 7 to split appears at rate \c\9 x ^ 2 ; a new cross or double bar that causes 7 and 
7' to merge appears at rate c\9~ x l 2 '^m- The rate for an existing cross or double bar to 

disappear is \c29 1 l 2 if 7 is split, and C2^ _1 ^ 2 |pT if 7 ancl l' are merged. Consequently, 
7 splits at rate 

Wl+^'^^lnL 2 (8.2) 

and 7, 7' are merged at rate 

LL 1 

(C 1 +C 2 )6»- 1 / 2 — EE r m iL'. (8.3) 

It is important that the constants c\ and C2 be the same for all loops and for both the split 
and merge events. This may seem a daring conjecture to make, but it has been verified 
numerically in [31 in a very similar model. It follows that the lengths of macroscopic loops 
satisfy an effective split-merge process, and the invariant distribution is Poisson-Dirichlet 
with parameter r s /r m — 9, see e.g. [BJ 29 and references therein. 

The case u & (0, 1) is different because loops split with only half the rate above. Indeed, 
the appearance of a new transition within the loop may just rearrange it: topologically, this 
is like ■<-» 8. This results in PD(|). Notice that this cannot happen when it = or when 
u = 1 on a bipartite graph. These considerations allow to formulate the following conjecture. 

Conjecture 2. Assume that d, 9, u, (5 are such that macroscopic loops are present. Then, 
as |A| — >■ co then k — > 00, 

• the random variable X^i=i f\X\ converges (in probability) to a constant, denoted v; 
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• the sequence of decreasing numbers (^p^y, ■ ■ ■ > ^p\A\ ) conver 9 es ft 71 probability) to a 
Poisson-Dirichlet distribution. More ■precisely, it converges to PD(0) if u = 1 and 
to PD(f) ifu€ (0,1). 

The case u = is a bit subtle as it depends on the graph. PD(#) is the right choice for 
bipartite lattices, while PD(|) should be expected otherwise. 

It turns out that Conjecture [2] is relevant for the discussion about symmetry breaking, 
even though the heuristics is rather indirect. If x and y are two vertices that are very far 
apart, the probability that they belong to the same loop is equal to the probability v 2 that 
they both belong to macroscopic loops, times the probability that they belong to the same 
element of the corresponding random partition. This can easily be calculated using Beta(i?) 
i.i.d. random variables Xi, X%, . . . , so that 



X u (1 - X,)X 2 , (1 - X 1 )(l - X 2 )X 3 , . 
has GEM($) distribution, which is closely related to PD(i?). Then 

P(£x, y ,o) « v 2 I>((1 - *i) 2 •••(!- X k ^fXl) 



v 2 



k>l 

$ 2 



EG 



, 4 + 2) (0 + l)(0 + 2) 



(8.4) 



_ V 2 

~ 0+1' 

The approximation should become exact in the limits |A| — >■ oo then ||x — y\ — > oo. 

8.2. Nature of pure Gibbs states. Let us focus on the case u = 1. The Hamiltonian 
Hj^ is (minus) the sum of transposition operators. Ferromagnetic product states of the 
form ® x |a), with a £ {—S, . . . , S}, are ground states: They are eigenstates of each T xy with 
eigenvalue 1; and this is the largest eigenvalue since T 2 y = Id. It is tempting to conclude 
that a ferromagnetic phase transition takes place (for d > 3) and that the pure Gibbs states 
are of the form (-)^ with £ S 2 : 

(■)n = lim lim (•)„•<!) , fc v * <? • (8.5) 
il h->0+ |A|->oo H A + h l^* u -b* 

The pure state (-)s 3 is represented by space-time spin configurations where long loops have 
spin 5, while finite loops have any spin values. It follows that 

(^)e 3 = vS, (Sl)s 3 = (Sl)s 3 = 0. (8.6) 

Decomposing the rotation-invariant Gibbs state into pure states, and using asymptotic clus- 
tering of pure states, we have 

(S x S y ) = ^(S x ■ S y ) — / (S x ■ Sy)fidQ 

Js 2 

3 (8.7) 

= $(Sx • S y )g 3 ps | ^^(5* )s 3 (S y )e 3 — \v 2 S 2 . 



On the other hand, using Theorem |3.3| and Eq. ( |8.4[ ) with $ = 2S + 1, we have 

(S^) = |5(5 + l)P(£ x ,»,o) « \v 2 S. (8.8) 
Eqs (8.7) and (8.8) agree in the case S = h. This should be expected, as is then the 



Hamiltonian of 



the Heisenberg ferromagnet (see Section 7.1). But the equations disagree 
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for all other values of S, in particular S = 1. Eq. (8.8) seems trustworthy as it relies on 



Conjecture [2j This suggests that the nature of symmetry breaking and the structure of pure 
Gibbs states are more subtle. Hopefully more light will be shed on these questions in the 
future. 

The case u = is similar, with the staggered magnetization replacing the magnetization. 
The calculations above confirm the existence of antiferromagnetic pure states when S = | , 
while the situation for S > 1 is less clear. In the case u 6 (0, 1) and S = A, we can check 
that Conjecture [2] is compatible with the breaking of the U(l) symmetry: Indeed, let S 1 
denote the unit circle in the plane 1-3; then 

^3 o3\ 1 / cl <~fl 1 c3 c3\ 



(S^Sy) — J (S^Sy + S^Sy) 



2 2lT / (S^Sy + S^.Sy}^dQ — J (S X S y ^ ^ S y ) 

Js 1 



(8.9) 



2 (S x )g 3 (Sy)g 3 + 2 (S x )e 3 (Sy)e3 



This is compatible with Theorem 3.3 (5|5*> = |P(^, y , ), and Eq. (8.4) with = 1 



Handwritten notes 

[Nl] J. Frohlich, Phase transitions and continuous symmetry breaking, Vienna lectures, available at 
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[N2] B. Toth, Reflection positivity, infrared bounds, continuous symmetry breaking, Prague lectures, avail- 
able at http://www.math.bme.hu/~balint/prague_96/ (1996) 

[N3] D. Ueltschi, Phase transitions in classical and quantum Heisenberg models, Tubingen lectures, available 
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